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Abstract. Motivated by asking how to combine parton showers with nonleading QCD 
matrix elements, we discuss a subtractive technique based on gauge-invariant Wilson- 
line operators and how this can be used to treat the soft region. 
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For precision phenomenology at future high energy colliders, it will be valuable to 
construct Monte Carlo event generators in which next-to-leading-order (NLO) QCD 
corrections can be incorporated in parton shower algorithms. There is at present 
no systematic method for doing this. Procedures based on subtractive methods 
^ ' have recently been proposed [1-3]. An important step in the implementation of 

this program is to show how to decompose Feynman graphs into sums of terms 
over different regions, with the terms arranged so as to correspond to factors in a 
factorization formula suitable for the Monte Carlo application. 

An example of such a decomposition is given in Ref. [2] for the photon-gluon 
fusion process in leptoproduction. In this case soft gluons do not enter at the 
leading power, so that leading regions do not overlap. To handle general cases, 
however, one needs to treat graphs with soft gluons and hence with overlapping 
leading regions. 

Fully numerical [4] or semi-analytical [5] subtraction methods have been devised 
to calculate NLO quantities that are infrared safe. These methods are not directly 
applicable in event generators that simulate the fully exclusive structure of the 
hadronic final states, since here the quantities being computed are not infrared safe 
in perturbation theory. In particular, one cannot use a cancellation of soft gluon 
contributions between real and virtual graphs. 

The technique we discuss in this talk is constructed so that the following prop- 
erties are satisfied [3]: 
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(a) The integrand for the hard scattering coefficient is to be an integrable function, 
even when the corrections are applied to a process that is not infrared safe. 

(b) The terms in the expansion of each Feynman graph should arise from matrix 
elements of gauge-invariant operators. 

(c) In particular, the necessary cut-offs on rapidity integrations should be defined 
gauge-invariantly. This involves the use of Wilson lines along non-lightlike 
directions [6]. 

(d) The evolution equations [7] with respect to these cut-offs should be simple, in 
the sense that there should be no power-law remainder terms. 

In what follows we illustrate this technique using a simple example, one-gluon 
emission graphs in leptoproduction. 

We consider a generic observable associated with the reaction j*(q) + q(p) — > 
g(k) + q(k'). We denote this by a[ip\, where ip is a weight function that contains 
the definition of the particular observable under consideration as a function of the 
final states as well as nonperturbative parts of the cross section, including the 
parton density. We denote by a and (3 the gluon's fractional lightcone momenta, 
a = k + /(xP + ), 13 = k~/(Q 2 /(2xP + )) (with x the Bjorken variable, Q 2 the photon 
virtuality, P + the hadron's plus momentum), and by the gluon's azimuthal angle. 
We represent a[<p\ as 
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a[<p]= da dp — - J(x, a, (3) M(a, (3) (p(x, Q 2 , a, (3, 0) . (1) 
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Here M is the next-to-leading-order matrix element. For the 7*g amplitude ob- 
tained by contracting the photon Lorentz indices with the projector corresponding 
to the structure function F 2 , M is of the form 
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FIGURE 1. The phase space of Eq. (1) in the a, f3 plane. 



(2) 



The factor J in Eq. (1) is a Jacobian factor [3] and defines the physical region in 
the a, (3 plane as the interior of the triangle in Fig. I. 

The infrared sensitive regions contributing to the leading power behavior of a[<p] 
are located in Fig. 1 as follows: the region in which the gluon is collinear to the 
initial state is a neighborhood of the axis (3 — 0; the region in which the gluon is 
collinear to the final state is a neighborhood of the axis a = 0; the soft region is a 
neighborhood of the origin a — 0, (3 — 0. The truly hard region lies away from the 
a = and (3 = axes. 

To obtain a decomposition for a into a sum of terms over the leading regions 
we employ the method of Ref. [6]. This generalizes the R-operation techniques of 
renormalization. See Ref. [8] for a related approach. We proceed from "smaller" 
to "larger" regions. For each region R, we remove the contribution from smaller 
regions, and construct an approximation to the matrix element valid in R, up to 
power suppressed corrections. Then we subtract any divergences that appear in this 
expression coming from larger regions. To ensure that the splitting between the 
terms is defined gauge-invariantly, at each step we demand that the counterterms 
be constructed from matrix elements of Wilson line operators, 



Vi(n) = V exp [ig J dyn-A(yn)^ , Vp{n) = Pexp (ig J dyn-A(y 
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with suitable directions n for the lines. Evolution equations in n enable one 
to connect the results corresponding to different directions. We introduce light- 
like directions p = (1,0, 0^), p' = (0,1, Or), and non-lightlike directions u = 
(u + ,u~,0t), u' = (u' + ,u'~,0t) (all with positive components). We define 77 = 
{2x 2 P +2 /Q 2 )u-/u + , and rj' = (Q 2 /2x 2 P +2 )u' + /u'' . 
The result for the term to be associated with the soft region reads [3] 
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Ms ^^ = ~p~ (a + rj'P) (3 ~ a {(3 + rja) ' (4) 

The first piece in the right hand side is just the soft approximation (a, (3 — > 0) to 
Eq. (2). The other two pieces are the subtractions associated with the collinear 
regions where a/ (3 or (3 /a go to zero. The first piece can be thought of as the one- 
loop contribution to the square of a vacuum-to-gluon matrix element of a product 
of eikonal Wilson lines taken along lightlike directions p, p' [6]. The subtraction 
pieces can be derived from operators analogous to those for the first piece, except for 
replacing one of the lightlike eikonal lines by a line along a non-lightlike direction. 
The parameters i], rj' act as gauge-invariant cut-off parameters. For example, in the 
second piece, which subtracts the divergence from the region (3 /a — > 0, 7/ provides 
a cut-off on the region of small a. 



By applying an analogous treatment we construct terms Mj, M F for the initial- 
state and final-state collinear regions. The term to be associated with the hard 
region is obtained by subtracting from the original matrix element the soft and 
collinear terms: 

M H (a, (3) = M-M S -M I -M F . (5) 

The explicit expression for M H is given in [3] . This expression is finite in all of the 
infrared regions, and independent of the choice of the non-lightlike directions u, u'. 

In the approach of [2,3], Mh gives precisely the subtracted hard-scattering func- 
tion to be used in the Monte Carlo to incorporate NLO corrections. The collinear 
terms correspond to factors in the showering. The soft term can be eliminated by 
a suitable choice of the directions for the Wilson lines: M s = for rjrj' = 1. Then 
all of the infrared contributions are associated with either initial-state or final-state 
branchings, as is appropriate to match the structure of the parton shower. 

It is interesting to consider the form of the collinear factors corresponding to 
the Monte Carlo scheme for the non-lightlike directions. In the symmetric case 
rj = rj = 1, the result for the initial-state collinear factor is [3] 

M rw)4 1+ ' 1+ f -i* (6) 

1 v ' f3 a (1 + a) a a + f3 w 

The first term in the right hand side corresponds to the standard form of the 
showering, with the coefficient of 1/(3 given by the quark — > quark splitting kernel. 
The second term, on the other hand, is non-standard. However, it is suppressed 
in the collinear region f3 — > at fixed a. That is, this form of the showering 
coincides with the standard one in the collinear limit and differs from it away from 
the collinear limit. The role of the second term is to provide an effective cut-off 
when a — > 0, i.e., a cut-off on the soft singularity of the splitting kernel. 

Observe that if one regulated the a — > behavior of the first term in Eq. (6) 
by subtracting its a — > limit, given by 2/(a(3), this would bring about an extra 
(3 — > singularity. This would not be suited for our application in a Monte Carlo 
algorithm. In contrast, the second term in Eq. (6) represents precisely what is, 
from our point of view, a better choice of a counterterm: it subtracts the a — > 
singularity without introducing any extra singular behavior at (3 — > 0. 

Finally, note that this counterterm cuts off the integration over the region of 
small a at a value of order (3. The procedure based on gauge-invariant subtractions 
tells us precisely where the cut-off is to be placed: the position of the cut-off on 
a turns out to be /9-dependent. In more physical terms, this indicates that the 
cut-off to be applied in the initial-state shower and the cut-off to be applied in the 
final-state shower are not to be set independently, but they are related. 

Acknowledgments. Thanks to J. Collins for collaboration on this topic. 
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